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Abstract 

We present exact combinatorial versions of bosonization identities, which equate the product of two 
Ising correlators with a free field (bosonic) correlator. The role of the discrete free field is played by the 
height function of an associated bipartite dimer model. Some applications to the asymptotic analysis of 
Ising correlators are discussed. 

1 Introduction 

The Ising model is a basic model of ferromagnetism in statistical mechanics. It has been extensively studied 
in two dimensions, at and around its critical temperature ([H]), where several powerful techniques are 
available: transfer matrices, pfaffian (dimer) representation, integrable systems and more recently discrete 
complex analysis and Schramm-Loewner Evolutions ([IS]). 

In the Ising model, the basic (order) observables are the spin variables; Kramers- Wannier duality (|37[l38]) 
maps these to disorder variables (|27j). In general, an order variable is a local random variable, while a 
disorder variable represents a local modification of the state space (eg [IB])- A classical object of study 
are correlators involving some order and/or disorder variables, especially in appropriate asymptotic regimes 

(PULED)- 

At the field theoretic level, the notion of bosonization introduced in and much developed afterwards 
(see in particular [13] and Chapter 12 in [T2]) expresses squares of Ising correlators as bosonic (free held) 
correlators. For the free field, the relevant order variables are electric insertions, and the dual disorder 
variables are magnetic insertions (eg [12 [Ml HE])- 

It has long been known that a planar Ising configuration may be represented by a dimer configuration on 
a related decorated graph ([2H[2H])j which may be analyzed through determinantal or Pfaffian techniques 
( 28J); these decorated graphs are not bipartite. Dimers on bipartite graphs are associated to a discrete 
height function, the asymptotic fluctuations of which have been extensively studied from the early 90's (see 
the survey [33]). In the type of scaling regime relevant here, these asymptotic fluctuations are described by 
a free field in a rather precise fashion. 

In the present article, we use (known) mappings and dualities between doubled Ising models, 8-vertex, 
6- vertex and dimer models (along the lines of [2 ), tracking order and disorder variables along the way, in 
order to establish exact combinatorial versions of the field-theoretic bosonization identities. These identities 
involve a pair of independent Ising configurations and a single bipartite dimer configuration. In combination 
with recent progress on the fine limiting behavior of dimer height fields f|17j). this enables to obtain the 
asymptotics of critical Ising correlators in the plane. Critical correlators in finite domains require some 
additional care and are the object of the upcoming [SI I19j. 
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The article is organized as follows. Bosonization identities are phrased in Sect. 2, along with a discussion 
of boundary conditions. Some consequences for the asymptotic analysis of critical Ising correlators are listed 
in Sect. 3. Relations with other approaches are described in Sect. 4. 

2 Bosonization rules 
2.1 Mappings 

Consider a graph T = (V,E) embedded on a torus £ = C/A, A = Z + rZ, 9r > (the planar case will be 
discussed afterwards). Let F denote the set of its faces and = (V\E^) denote its dual graph, so that 
V* ~ F, ~ E (if e £ £, e f £ denotes its dual edge). To each edge e g E we associate a coupling 
constant J e > 0. A configuration of the Ising model consists in an assignment of a spin a v = ±1 to each 
vertex w of L. By planar duality, one can equivalently assign spins to faces. The weight of a configuration 
(<?v)vev i s: 

w((a v ) veV ) = exp(-2/3 Jgl^^J 

e=(W)£-E 

where the inverse temperature f3 is a fixed positive constant. The weight is invariant under global spin flip 
(<7„)„ <-> (— <7„)„. In the low temperature expansion, one represents a configuration (cr„) by the even degree 
subgraph (polygon) P„ = (V\El) of rt = (V\E*), where e = (W)t € £ ff iff cr„ ^ ov- Clearly 

= n w ( e ) 

if we set 

w(e) = exp(-2/?J et ) g (0,1] 

A polygon P has even degree at each vertex. Given P, it is easy to see that the parity of the number of 
edges crossed by a closed cycle on T depends only on the homology class of this cycle (in i?i(£,Z)). Thus 
to any polygon one may associate signs (ea{P),£b{P)) € {±1} 2 where ea = +1 or —1 according to whether 
an A-cycle crosses an even or odd number of edges of P, and similarly for ub- Plainly, a polygon P comes 
from a (periodic) spin configuration a iff (ea{P), Eb(P)) = (1, 1). (We could frame this discussion in terms 
of discrete 1-forms with values in {±1}). 

At this point it is natural to introduce spin configurations with periodic (p) or antiperiodic (a) boundary 
conditions corresponding to A and B periods. One way is to consider spin configurations a on the lift of T to 
C/(2A) such that cr(x + 1) = e a o-{x), a{x + t) = e b o(x), where (e A ,e B ) = (1, 1), (-1, 1), (1, -1), (-1,-1) 
corresponds to (pp), (ap), (pa), (aa) boundary conditions. Associated to such a spin configuration a, we have 
a polygon P a , where now (ea, Eb)(P<j) may take any value in {±1} 2 . 

Thus consider the space of polygons on V with weights 

w(P) = ] [ w(e) 

eEEp 

The state space is partitioned in four blocks corresponding to (pp),. ■ ■ ,(aa) boundary condition. Each polygon 
lifts to two spin configurations ±er, in a measure preserving way (up to normalization). 

Let us also introduce the medial graph M = (Vm, Em) of T (or derived graph, in the restricted context of 
planar graphs): the vertices Vm of M are set at the midpoints of edges of T; two vertices of M are adjacent 
if they correspond to consecutive (in cyclic order) edges around a vertex of L. Note that M is 4- regular (all 
vertices have degree 4) and may be identified with the medial graph of 1^ . Every edge e of M corresponds 
to two edges of L that have a common endpoint v £ V and lie on the boundary of a common face / g F. 
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Figure 1: (portion of a) graph T; dual graph T^; medial graph M; bipartite graph C 



This establishes a natural correspondence between edges of M and pairs (v, /) € V x F where v is a vertex 
on the boundary of F. 

Given a 4-regular planar graph (such as M), we define an 8-vertex configuration (or 8V configuration for 
short) as an orientation of the graph such that the number of incoming (or outgoing) edges at each vertex 
is even. Around each vertex, there are 8 possible configurations (orientations of the four adjacent edges). 

From dual Ising models to the 8V model. 

Let us now consider two spin configurations a and a\ with a defined on V and er^ defined on F ~ 
with the same boundary condition in {(pp), . . . , (aa)}. To the pair (a, er^) we associate an 8V configuration 
on M as follows: Let e be an edge of M corresponding to the pair (v,f) € V x F. If a(v) = cr'(f) (resp. 
o~(v) — —cr'(f)), we orient e such that v is on its right handside (resp. left handside). In other words we have 
a reference orientation of M where vertices of T are on the righthand side of oriented edges of M. Other 
orientations are associated with an edge spin configuration: v(e) = a(v)cr'(f), where e € Em corresponds to 
(v, f) <E V x F (see Figure [2]). The orientation (8V configuration) associated to v agrees with the reference 
orientation on e iff v(e) = 1. 

Let us now assume that (a(v)) v ^v and (f (/))/eF are sampled independently from Ising distributions 
(with couplings {J e )eeE, (^ e t) e teBt an d associated edge weights (w(e)) e , (w(e>)) e -t) and same boundary 
condition in {(pp), ■ ■ ■ , (aa)}. The weight of the configuration (a, <j^) (up to normalization) may be expressed 
as a product of factors w(e)w(e^), where e, is a pair of dual edges of T, that corresponds to a vertex of 
M. The 8V weights at this vertex are as follows (w, is the weight of the i-th vertex configuration, see Figure 
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Figure 2: Graph T (solid) and (dashed); medial graph; spin variables and associated 8V configuration 
type: 



1 


3 


5 


7 


2 


4 


6 


8 



uj\ = lo 2 = w(e) = ujq = 1 

ll>3 = uii — w(eJ) 1M7 — lc>8 — w(e)w(e < ) 

Remark that (a, er^) and (— a, — a^) project to the same 8V configuration, while (a, — cp), (—a, a 1 * ) projects 
to the configuration with all arrows reverted (it has the same weight). When 8V weights are invariant under 
reversal of all arrows, we denote u>i2 = u>i = u> 2 etc. 

From the 8V model to the 6V model. 

These weights may be written as a function of the medial edge spins v{e m ) — a(v)a' (/). If we number 



14 n 

2 g the medial edges in Figure 21 we may write the weight as 



u>(ui, ...,U4,) = a(l + V1V2V3V4) + b{v\v 2 + + c(y\V± + v- 2 v- i ) + d{v\v z + v 2 v A ) 

because of the invariances ■ ■ ■ , V4) = u>(—vx, . . . , —^4) and cd(i>i, v 2 , ^3, V4) = ^(^3, V4, Vi, v 2 \ Note that 
ijj{v\ 1 . . . , 2/4) = if V\ . . . V4 = — 1, which enforces the 8V condition. Solving 

2(a-b + c-d) = ui(l, -1, -1, 1) = W12 

2(a + b-c-d) = 1, -1,-1) = w 34 

2(a + b + c + d) = cj(1, 1, 1, 1) = w 56 

2(a-b-c+d) = w(l, -1,1,-1) = w 78 

we get 

8(a, 6, c, d) = (W12 + w 34 + oj 56 + w 78 , -W12 + ^34 + w 56 - w 78 , w i2 - a; 34 + w 56 - w 78 , -u) 12 - w 34 + w 56 + w 78 ) 



The abelian duality for the 8-vertex model is obtained as follows. Let us start from the 8V partition 
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function 

z = e n^M 

where LO e {v) is the 8V weight at e, which is written as a sum of characters of {±1} 4 - Expanding Ile^et^) 
in monomials in the v e variables, we notice that only monomials with even degree (0 or 2) in each variable 
contribute to the partition function. To each contributing monomial, we associate a medial spin configuration 
(i) e )eeB M as follows: i> e — \\i the partial degree of v e is and z> e = —1 if the partial degree of v e is 2. This 
yields 

Z = Z 

where Z is the 8V partition function with weights 

U>12 — W34 + W56 — LO-js 

4 

— Wi2 + W34 + W56 — W78 
4 

U!\2 + W34 + W56 + W78 
4 

— 0>i2 — W34 + W56 + W78 
1 

Note that the weight mapping (W12,...) h-> (cDi2,..-) is involutive. The spin variables i/ e , e € -Em, arc 
natural order variables. It is easy to see that under duality, they are exchanged with disorder variables 
defined as follows. A disorder variable £ e at the medial edge e <G Em is a defect splitting it in two half-edges 
with opposite orientations. It is thus a modification of the state space, which we still denote as a random 
variable by (a standard) abuse of terminology. The duality identity then reads: 

(z/(ei) . . . ^(e m )£(e m+ i) . . . £(e m+ „)) w = (£(ei) . . . £(e m )i/(e m+1 ) ...v(e 

m+n ) ) w 

where the LHS is Yl u i e i) ■ ■ ■ v { e m) ELeE w eWi where the sum bears on configurations with defects at 
e m+ i, . . . , e m+n . The RHS is defined symmetrically. 

In presence of a disorder £ e , there are two opposite spin variables v{e~), v{e + ) corresponding to the two 
half-edges e + ,e~ of e. Then duality maps the pair i/(e + )£(e) to £(e)f(e~), and vice-versa. 

Assume that the two sets of Ising weights {w{e)) eeE , (w(et)) e t e st satisfy the Kramers- Wannier duality 
relation: 

10(e) + w(e f ) + w(e)w(e f ) = 1 (2.1) 
for each pair (e, e^) of dual edges. Then the associated 8V weights are 

(W12, W34, w 56 , w 78 ) = (10, w', 1, W) (2.2) 

with w = w(e), w' — w(e') for short. In any 8V configuration on a toroidal graph, the number of sinks (type 
7) equals the number of sources (type 8) . Consequently, we get the same configuration weights if we change 
the local weights to: 

(W12, W34, w 56 , uj 78 ) = (w, to', 1, -ww') (2.3) 
Applying duality, we get the same partition function with the dual weights: 

(Wl2, CJ34, W 56 , (2)7$) — ^(w — w 1 + 1 + ww', — w + w' + 1 + ww' , w + w' + 1 — ww 1 , — w — w' + 1 — ww') 

= ^(1 — io', 1 — w, 1 — ww', 0) 

This set of weights defines a 6- vertex model (as sources and sinks get zero weight), with weights: 

(a>i2,W34,a> 56 ) = -(1 - w', 1 - w, 1 - ww') (2.4) 



&12 = c/2 = 
W34 = 6/2 = 
&56 = a/2 = 
w 78 = d/2 = 
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Moreover, 

LJl2 2w W34 1 — W 2 

a) 56 1 + w 2 ' w 56 1 + w 2 
so that cD| 6 = d)j 2 + 0)34, ie these are the weights of a 6V model at the "free fermion" point. 
Prom 6V to bipartite dimers. 

Fan and Wu showed ([20]) that the 8V model on the free fermion line, ie satisfying the vertex weight 
relation: 

U1LO2 + ^3^4 = ^5^6 + ^7^8 

can be mapped exactly to a dimer model on a decorated graph, in a way rather similar to the Temperley- 
Fisher mapping of the Ising model to a dimer model. In the case of the 6V model at the free fermion point, 
one can find a bipartite dimer representation, which we now describe. 

We start from a 4-regular planar graph M (plainly, this works for toroidal graphs), with checkerboard 
coloring of faces. A decorated graph C is obtained as follows. Each vertex v of M is replaced with a city (as 
in the urban renewal metaphor, see |35j). a quadrangle with a vertex on each edge abutting v. Each city has 
four internal edges and two adjacent cities are connected by a road edge. Formally, vertices of C correspond 
to pairs (v, e), v £ M, e e Em abutting v; ((v, e), (v 1 , e')) is an edge of C iff either e = e' = (vv') (road) or 
v = v' and e, e' are two consecutive (in cyclic order) edges around v (city street). 

Plainly, C is bipartite: the city faces are quadrangles, and the degree of other faces is twice the degree of 
the corresponding face in M (an additional argument is needed to ensure consistency when wrapping around 
the torus; this is obtained from the checkerboard coloring of the faces of M). A dimer configuration or perfect 
matching on C is a subset m of edges of C such that each vertex of C is the endpoint of exactly one edge 
in C. Given a set (w e ) e< =E c °f positive weights associated to edges of C, one defines a dimer configuration 
weight by 

w ( m ) = H w e 

For background on the dimer model, see eg [53]. 

We can choose the checkerboard coloring of faces of M and the two-coloring of vertices of C in such a 
way that if (bw) is a road edge of C oriented from black to white, the face of M on the RHS (resp. LHS) is 
black (resp. white). A 6V configuration consists on a orientation of edges of M (with the 6V rule enforced 
at each vertex). Given such a configuration, each road edge of C inherits an orientation. We decide that 
(bw) £ Eq belongs to the corresponding dimer configuration iff (bw) is oriented from black to white. This 
completely determines the perfect matching, up to a local ambiguity for each type 6 city, where two opposite 
city streets are matched. See Figure [3] for a graphical representation of the correspondence. 

It is obvious that the correspondence is measure-preserving if the 6V weights are given by 

(wi2,w 34 , uj 56 ) = (sin 9, cos 6, 1) 

(where 9 may depend on the city; the weight of the corresponding edge is w = tan(0/2)) and the dimer weights 
are 1 for each road, sm9 for horizontal city streets and cos# for vertical city streets (in the coordinates of 
Figure [3]). The various graphs derived from T are illustrated in Figure [Tim the generic case and in Figure [4] 
when r is the triangular lattice. 

Associated to the 6V model or its bipartite dimer representation is a height function <f> which may be 
described as follows. To each face / of M one associates a height </>(/) such that 4>(f) — 4>(f) — 7T (resp. 
— 7r) is crossed from left to right by an oriented edge in the 6V configuration. The 6V rule ensures that 
5Zi=o(/i+i — fi) — 0, where /o, • • • , /3 are the cyclically ordered faces around a vertex v £ M. This gives a 
consistent definition of (up to an additive constant) on a fundamental domain. On the torus, cj) is additively 
multivalued. 
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Figure 3: local 6V configurations on a 4-regular graph (type 1,2;3,4;5;6) and associated local dimer configu- 
rations on a decorated graph (bold edge: dimer) 

In terms of the height function, there are two natural (at least in the present context) types of insertions: 
electric and magnetic charges. If / is a face of AI, the r.v. e %a ^> is referred to as an electric charge a 
inserted at /. A pair of opposite charges e la ^ )e~ la ^^ does not depend on the choice of a base value for <j>. 
We may choose the additive constant of <fi so that (/>(/) = mod 2ir if / is a black face of M (corresponding 
to a vertex of T); then <p(f) — tt mod 2ir if / is a white face (corresponding to a vertex of T'). If / G T, 

e |0(/) =e -m) =cos(0(/)/2) 

and if / G rt, 

-ieiHf) = ie-iW) = sin(0(/)/2) 
An order line is a product Ileei ^( e ^' wnere 7 is a path from / to /' on M*. We have: 

°UW) = II = e 4(0(/,) " 0(/)) 

provided /, /' are of the same color (in the bipartite coloring of Aft). Finally we get the identification: 
a(f) = cob(0(/)/2) if / G r and a{f) = sin(0(/)/2) if / G rt. 

Inserting a disorder variable at e G E^j amounts to introducing a magnetic charge ±1: the height function 
becomes locally multivalued, picking an additive constant ±7r (depending on the type of defect: source or 
sink) when cycling counterclockwise around the defect. Remark that this has always to be compensated by 
another (opposite) magnetic charge elsewhere. 

In terms of the dimer mapping, disorder variables may be interpreted as monomers ([23 ). In the local 
correspondence, a 6V configuration with a sink defect (magnetic charge -1) at e G Em — (bw) G Ec 
corresponds to a modified dimer configuration in which b is matched to e (the midpoint of (bw)) while w is 
matched in its city. Conversely, a source defect corresponds to w being matched to e and b matched in its 
city. Removing the half-dimer, one gets the more classical formulation of a monomer defect at b or w . 

Specifically, if (.)dimer denotes the unnormalized expectation (ie (l)dimer is the partition function) under 
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Figure 4: Triangular lattice T (solid) and its dual hexagonal lattice (dashed); associated 4-regular graph M 
(Kagome lattice); associated bipartite graph C (4-6-12 archimedean lattice) 

the dimer measure on C, we denote 



the unnormalized dimer measure on C \ {&i, tui, . . . , b n , »)„}, where bi, . . . , b n (resp. wi,..., w n ) are n black 
(resp. white) vertices on C (using a magnetic operator notation, where as usual disorder variables represent 
a modification of the state space). If the number of black and white monomers does not match, we set 
riiLi ^M&») IlJLi ^-iC^j) — ( as the state space is then empty). 

With these conventions, we have £(e) = O x {b) + 0-i{w) and v{e~)(,{e) = -v(e + )£(e) = O x {b) - Oi(w). 
where e~ is the half-edge of e containing b. 

Correlators. 

We are now concerned with tracking correlators of order and disorder variables through the mappings 
(Ising, dual Ising)^ 8V — > 8V^ 6V^ dimers. 




dimer 
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In the Ising model on T = (V,E), the basic order variable is a spin variable a v . By Kramers- Wannicr 
duality, such order variables correspond to disorder "variables" in the dual Ising model. Following Kadanoff 
and Ceva, a pair of disorder variables fiffif is represented by a simple path 7 from / to /' on (a disorder 
line); the weight of a configuration is then modified to 

w'{(a v ) veV ) = Yl w ^ IT w ( e ) 

A general order-disorder correlator for the Ising model on T is written as: 

/ 2n 2m \ hv 

(iI^IlM/i)) 

\i=l j=l I r 

where hv E {p, a} 2 designates the periodicity type of the spin configuration, and the /j's are connected 
pairwise by disjoint disorder lines drawn in a fixed fundamental domain (associated to the choice A, B of 
homology basis). We also assume that all insertions i;,*, fj are disjoint (at no cost since a 2 = 1 and /u? = 1). 

Let us consider simultaneously a dual Ising model on T\ with a general correlator of the form: 

I In' 2m' \ h v 

(n^n^K)) 

\i=i j=i I rt 

For definiteness, let us assume that the spin sites Uj (resp. /•) are paired by simple "order lines" on V (resp. 
r^), and that for each of the two configurations, order and disorder are disjoint (and are also disjoint from 
A, B cycles bounding a fundamental domain). 

It is clear that in the mapping (Ising, dual Ising)— > 8V, an order line from, say, v\ to vi on T can be 
represented by an order line, ie a simple path 7' from v\ to «2 on M* . This is simply saying that 

a(vx)a(v 2 ) = Yl K et ) 

Let us now consider an oriented disorder line 7 from f\ to f% on T\ say. One way to represent it goes 
as follows: let e be an oriented edge on 7, which separates the site v on the left from the site v' on the 
right. We add a new site (ye) on e\ close to v, and set a spin variable a(ve) = —a(v). Then the factor of 
the configuration weight corresponding to (vv r ) is e P J e cr (. ve )< J ( v ) (rather than e PJe< y ( v )< J ( v ) m the absence of 
disorder). 

Correspondingly, if e is an edge on M which ends on 7 and is on its LHS (corresponding to v in V and 
/ in F), we define two opposite edge spin variables: a(v)a^(f) and <j(ve)a'(f); in other words, we insert a 
disorder variable £(e). Consequently, we can represent the Ising disorder line 7 on rt as an 8V disorder line 
7 on M\ which tracks 7 on its LHS. See Figure [5] for an illustration of the situation. 

Plainly, we can recover the spin configuration on M* from the 8V configuration up to global spin flip. 
Thus we denote by <r(/) (resp. /i(/)) the endpoint of an 8V order (resp. disorder) line, where / is a face of 
M, corresponding to a vertex of r or T\ For given endpoints, changing the pairings of insertions or moving 
the disorder lines may only change the sign of correlators. 

Hence we obtain the first identity: 

/ 2n 2m \ hv 1 2ri 2m' \ I n+n' m+m' \ 

(n^n <"(/.;•)) (n^n^K-)) = 2 w^nn^ii n^) 

\i=l j = l I p \i=l j = l I p t \ i=l eg7i j = l eeY, I „., 

I 2n 2n 2 m 2m \ 



2 ( i eA=KEB = v n «r(«o n ^(/i) n ma n *w 



1 8V 

(2.5) 
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Figure 5: Left panel: square lattice (solid), primal and dual spin configurations, disorder line (dashed, 
oriented from bottom left to top right). Right panel: corresponding 8V configuration with disorder line 
(dashed) 



Here the 8V weights are as in (2.2); the 7,'s (resp. 7J) are order (resp. disorder) lines on M* pairing the 
order (resp. disorder) insertions v\, V2 n , /(> fyn' ( res P- /i> • • • > hm-, v 'ii ■ ■ ■ i v 2m')- The signs e^es are the 
products of 8V edge variables v along A, B cycles on M'. We identify p = 1, a = — 1 for periodic and 
antiperiodic boundary conditions. The factor 2 accounts for the spin flip symmetry (a, & 1 ) <-> (— a, — a*). 



The next step is to change the 8V weights from (2.2) to (2.3). This is somewhat complicated by the 
presence of disorder lines. One can decompose an 8V configuration with disorders (as in Figure [5]) as a 
collection of oriented curves: four curves are created at a source (type 8), four are destroyed at a sink (type 
7); at a disorder insertion, two curves are created or two at destroyed. Denoting by p + , p_ the number of 
edge sources and sinks, N% the number of vertices of type i, we have: 

4N S + 2p + = 4N 7 + 2p_ 

Note that p + = p- mod 2, since all our disorder lines on (ie those coming from Ising disorder lines on 
r or r^) have even length. Hence 



where 8V weights on the LHS (resp. RHS) are given by (|2.3| (resp. (2.2)). It is easily checked that if 7 is a 



disorder line, P+(7),P-(7) the number of edge sources/sinks along this line, and 7 is an (edge disjoint) order 
line with the same endpoints, then: 

P-lCt)— P— (7) t— r , 

(-1) — a — =n^ f ) 

Consequently: 



/ 2n 2m \ 

\i=l 3=1 I 



1 2n 2m' \ V j n+n' m+m' 

n^^n^K)) =2(i eA =Ma=« n n n^n^) 

\i=l j = l I p t \ i=l e£7i 3=1 eGy^ ee^/j 

I 2n 2n' 2m 2m' \ 

= 2 ( i EA=h , EB=v n n n m/x^o n mk>ko 



8V> 



i = l 



3 = 1 



3 = 1 



(2.6) 



8V 



where jj is an order line with the same endpoints as 7', and the 8V weights are as in (2.3) 



10 



Now the 8V duality simply exchanges order and disorder insertions, so that: 



hv 



I n-\-n 



m+m 



/ 2n 2m \ ,bU / 2ri 2m 

e (n^)n^)) (n^^n^K)) =2(nn^nn^)ii^)) 

h,v=±l \i=l 2=1 / r \i=l 2=1 / pt \ i=1 e£ Ti 3=1 e 67l- e^Ti / 

/ 2n 2n' 2m 2m' \ 

= 2 ( n n n ^w-) n * w^) ) 

3=1 3=1 / 6 y 

(2.7) 



i=l i=l 



where the 6V weights are given by (2.4). 



2.2 Spin structures 



We would like to obtain a bosonization identity of type (2.7) which a factorized LHS. For this purpose, we 



need to discuss in greater details the relation between boundary conditions in the various models considered. 

First we carry out the discussion in the absence of insertions. To an 8V (a fortiori 6V) configuration, we 
can associate a pair of signs (ea,£b), where £[ 7 ] = Ilee7 ^( e ^) e where 7 is a closed cycle on and 

[7] is its homology class in H 1 (E,Z/2Z) ~ i?i(£, {±1}) ~ {1, -l} 2 . By the 8V condition, the RHS depends 
on 7 only through [7], and we have £[ 7l ] + [ 72 ] = £[ 7l j£[ 72 ]. 

We have 

(1>6V = (l) 8 v - 



E < 

h,v=±l 



1 \ hv / 1 \ hv 

l )r Wrt 



and would like to suitably twist this identity by a character of ffi(£,Z/2Z) (a spin structure). If 7,7 are 
two homotopic simple cycles on M\ we have: 



n^( et )n^ et ) 



jji/( e t) 



(%]) 



W/6V 



\e£7 



eS7 



8 V" 



\eg7 



8V 



\e67 



6 V 



For a spin configuration a (with periodicity conditions in {±1} 2 ~ {p, a} 2 ) on r or , we set £r 7 i (a) = £7] (P) 
where P is the corresponding low- temperature polygon; this depends solely on the boundary condition. This 
gives an interpretation of (closed) order cycles: £[ 7 ](cr) = \\. e= i vv i\^ cr(v') a(v) (in the presence of possibly 
antiperiodic boundary conditions, each factor in the RHS is unambiguously defined). If (h,v) G {±1} 
represents the periodicity conditions, ea — h and Eb = v. 

Dually, we can consider a disorder cycle 7 ( "antiferromagnetic seams"), say drawn on T: the coupling 
for the Ising configuration on is negated on each edge crossed by 7. By drawing 7 on the boundary 
of a fundamental domain, introducing a disorder cycle is equivalent to shifting the periodicity condition. 
More precisely, let us identify £^(£,{±1}) with {±1} 2 by setting [A] = (-1,1), [B] = (1,-1). We define 
an involution [7] i-> [7]* by [A]* = [B], [B]* = [A]; and an F 2 ~ {±l}-bilinear pairing by [A] A [B] = 
[B] A [A] — —1, [7] A [7] = 1 for all 7. Then inserting a disorder cycle 7 shifts the periodicity type from (h, v) 
to (h, v) + [7]*; and inserting an order cycle 7 introduces a sign [7]* A (h, v). 

Now let us consider a spin configuration on T with an order cycle 7 and a spin configuration on with 



a disorder cycle 7. As in (2.7) (except we have closed cycles instead of open paths), we have 



(%])6v=(n^ et )n^ et ) 



8V 



h,v=±l 



\ £ Wt Wrt 



E 

h,ti=±l 



'M* A (h,v)) (l)p" (l)p" + ' 7 '* 



ll 



Let q : {±1} 2 — > {±1} be the quadratic form given by q(l, 1) = 1, q(a, b) = — 1 otherwise. For (/x, v) £ {0, l} 2 
(a sector), we set: 

h,v=±l 

so that A " / (-)r is a signed measure with state space the disjoint union of the four spin configuration spaces 
corresponding to the periodicity types in {±1} 2 . Remark that (pairs of) order and disorder variables are 
defined consistently in these four spaces; and that in the low-temperature representation, there is a 2-1 
correspondence between this total configuration space and the space of polygons (even degree subgraphs of 
r). Besides, we have the inversion formula 

(■)r v = \ E q(h,v)hWri-)r) 



Applying the previous identity to [7] = A, B, A + B, we obtain: 

2(1) 6V = (00)(00)t + (10)(10)t + (01)(01)t + (ll)(ll) t 

2(e A ) 6V = (0Q)(01)t - (10)(ll)t + (01)(00)t - (ll)(10)t 

2(e B ) 6V = (00)(10)t + (10)(00)t - (01)(ll)t - (ll)(01)t 

2(e A e B ) 6V = (00)(ll)t - (10)(01)t - (01)(10)t + (ll)(00)t 

where for legibility we write (hv) (resp. (hv) r ) for (1)^ (rcsp. identifying {±1} with Z/2Z. Then a 

linear combination of these relations yields: 

2(g(£A,£ S )) 6 v = 00 (l)r 00 (l)rt 

and more generally 

2(q(e A ,e B )e A e B ) 6V = ^(l) r ^(l) r t 

for n,u G {0,1}. 

The argument may be carried out in the presence of insertions, as in ( |2.7[ ). Indeed, as long as order 
and disorder lines are drawn in a fundamental domain, and the spin structures/boundary conditions are 
materialized by cycles on the boundary of the said domain, we have (with the 7's as in (2.7|) 

^ v I 2m \ W I In' 2m' \ / n+n' rn+m' ' 

(n^oiiM/i)) (li^i^u^^n = 2 ^^nn^ii n^n^) 

\i— 1 j — 1 I p \i— 1 j — 1 I pf \ i—1 e^7i j — 1 e^7^ e ^7j ' 

/ In In' 2m 2m 

= 2 / q(s A , e B )e A e B J[ fx( Vi ) [] J[ a{fMfj) II < v 'M v 'i) 

\ i—1 i—1 j—1 i — 1 

(2.8) 

for (j,, v E {0,1}. 

The reader will have noticed the close parallel with spin structures for dimers |29j . In higher genus, this 
has been completely clarified in |10[ lllj . in particular in terms of the Arf invariant, which allows to extend 
the above discussion to the case g > 2. 



2.3 Boundary conditions 

In simply-connected domains, we avoid difficulties related to spin structures; however a discussion of bound- 
ary conditions is needed. 
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Standard boundary conditions for an Ising configuration are the wired and free conditions. A wired 
boundary arc is one on which spins are constant; this may be enforced by setting edge weights to zero along 
this arc. On a free arc, there is no interaction with the outside; for symmetry (wrt duality), this may be be 
represented by adding an outside vertex connected to each vertex on the free arc by an edge with weight 1. 
We will consider the case of a finite planar graph (thought of as the approximation of a macroscopic simply 
connected domain) with outer boundary partitioned into alternating wired and free boundary arcs (thought 
of as approximating finitely many components of the boundary). 

For simplicity, let us first consider the case where the boundary of the primal graph T is partitioned in 
two (possibly empty) arcs which are respectively wired and free, as illustrated in Figure [6] Correspondingly, 
the dual graph has a free and wired arc. 




Figure 6: Left: graph T with a wired arc (bold); dual graph (dashed) with a wired arc (bold). Right: 
corresponding 4-regular graph M 

Because the spin is constant on a wired edge, we can construct the 4-regular graph M (carrying corre- 
sponding the 8V model) in such a way that edges of M bounce off wired edges of T or F^ . Then each vertex 
of M corresponds to a pair of dual edges with nondegenerate weights. 

In the mapping to dimers, one still replaces each vertex of V with a city; and these cities are connected 
by a road or a chain of roads (with weight 1) along each edge of M. Edges of M now include one edge per 
wired edge on T or (properly counted if the boundary arc is not simple) , and one extended (macroscopic 
edge) for each pair of changes of boundary conditions. Remark that each change of boundary condition from 
wired to free or free to wired corresponds to a distinguished vertex of the dimer graph C, and that these 
distinguished vertices are the endpoints of the external road in C. This procedure produces a dimer graph 
C which may contain chains of degree 2 vertices (in faces of T or which have several consecutive edges on 
a wired boundary arc). The dimer model is essentially unchanged if such a chain is replaced with another 
chain with the same number of vertices modulo 2, see Figure [7J 




Figure 7: Left: a face of T adjacent to a wired arc (bold). Middle: corresponding graph C (locally). Right: 
an equivalent dimer graph. 
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The situation is somewhat more involved when the number of boundary arcs (again alternatively wired 
and free) is 2n > 2. Then we have the option of adding external wirings connecting some of the wired 
components; this external wiring materializes a given partition of the set of the n wired boundary components. 
For instance, for two wired and two free boundary components, one can consider two cases, depending on 
whether or no the two wired components are connected by an external wiring. 

In order to preserve planar duality, we consider the case where the external wiring can be realized by 
disjoint edges, or equivalently the partition of the wired components is non-crossing. As is well known, there 
is a duality correspondence between the non-crossing partitions of the n wired arcs and the non-crossing 
partitions of the dual wired arcs, and also a correspondence with the non-crossing pairings of the 2n endpoints 
of the boundary arcs. These different configurations are enumerated by Catalan's numbers. It is easy to 
check that a choice of non-crossing partition for the wired components on V corresponds to the choice of a 
non-crossing pairing of the 2n distinguished vertices on C (located at the boundary condition changes) by 
external extended edges. See [HI [36l |8] for related considerations. 

To summarize, given a (finite, connected) planar graph Y with an outer boundary consisting of 2n 
alternating wired and free arcs, and a non-crossing partition of the wired components, we can construct a 
4-regular graph M carrying a 6- vertex model and a bipartite graph C carrying a dimer configuration in such 
a way that: 

/ 2n 2m V 1 2n' 2m' V J n-\-n' m+m' \ 

(n^n^/i)) (n^n^K-)) = 2 iin^n u^u&n ^ 



= 1 e£7i j = l e g 7 ' ee^j 



6V 



where the 6V weights are given by 2.4 and order and disorder lines 7, 7', 7 are as in the toroidal case. This 
also applies when the whole outer boundary consists of a single wired or free arc (in which case there is no 
external edge in C). 

It is also standard to consider ± boundary conditions, where the spin values on some boundary arcs are 
set to +1 or — 1. This removes the spin flip invariance of the model, which is inconvenient for our purposes. 
Note however that these conditions may be realized using wired and free boundary conditions. For instance, 
if (•)_)_ (resp. (-)u>) represents the unnormalized expectation for the Ising model on a planar graph Y with + 
(resp. wired) boundary condition, we have 

l + a(b) 



where b is any boundary vertex of Y. Indeed, l,r(b)=i — 1+ ^ b ^ . Note that boundary edge weights on a ± (or 
wired) boundary component are not counted in the spin configuration weight. 

Similarly, if the boundary is split into a + and a — boundary arc, this may be represented in two ways: 
by considering a wired- free- wired-free boundary (and no external wiring between the two boundary arcs), 
with each free "arc" consisting of a single edge at the endpoint of the ± arcs, and then: 

i + g(fe+) l-g(fe-) 
2 ' 2 

where 6+ (resp. b—) is a vertex on the + (resp. — ) boundary component. Alternatively, one can consider 
the boundary as a single wired component and place a pair of disorder variables at the boundary condition 
changes. 



2.4 Bosonization rules 



In the identities (2. 8), (2. 9), we have a relation between a product of two Ising correlators involving order 



and disorder operators, thought of as local operators, and some 6V (or equivalently dimer) quantity, which 
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might involve macroscopic disorder lines. This seems rather unpractical, so we will focus on the case where 
there is no such macroscopic line (order lines correspond to electric correlators at the endpoints). Each spin 
variable (on the primal graph F or its dual graph F*) corresponds the endpoint of a disorder line on the 6V 
model; each (primal or dual) disorder variable corresponds to the endpoint of an order line and in disorder 
line in the 6V model. Consequently, in these identities, if the insertions (ie order or disorder variables on 
the primal or dual graph) Vi, fj, f[, v'j are grouped in microscopic blocks of even cardinality, then the RHS 
is expressed in terms of local 6V (order and disorder) variables. 

In the simplest (most interesting) case, such a microscopic "block" consists in v e T, / e adjacent. 
This corresponds to an edge e in the 8V/6V graph M and a road (bio) in the dimer graph C. We denote 
e~, e + the two half-edges of e, with e~ corresponding to b. We recall that 0\{b) (resp. 0-\{w)) represents 
a monomer insertion at b (resp. w), thought of as a magnetic insertion. We also denote by v, f the two 
corresponding faces of M and by <j> the 6V height field which is defined (at least as a multivalued function) 
on the faces of M. Finally, ip(vf) — a(v)fi(f) is the Kadanoff-Ceva fermion. This yields the bosonization 
rules recorded in Table [T] 



Table 1: Bosonization rules 



Ising x Ising^,8V 


8V 


6V 


Dimer 


<j(v)aHf) 


a(v)a(f) 




1 (b) + 0- 1 (w) 


m(/)m t («) 


»(v)a(vMf)a(f) 




±{O x {b)-0_ x {w)) 




H(v) 


a[v) 


cos(^(«)/2) 


* 1 (f)Kf) 


tif) 


*(/) 


sm(0(/)/2) 


i,(vf) 


a(v)a(f)»(f) 




(Oi(6) + O_iH)Bin(0(/)/2) 




cr (/) cr («)A*(^) 


H(f)n(v)<r(v) 


(0 1 (6) + 0_ 1 H)cos(^(u)/2) 




m(/)m(«) 


a(f)a(v) 


cos(0(«)/2) rin(0(/)/2) 



Here ±(Oi(b) — 0_i(w)) is taken to be +(O x (b) — 0_i(w)) at the start (say) of an 8V disorder line and 
(Oi(b) — (D-i(w)) at the end (this comes from the fact that the sign of £(e)^(e ) depends on the half-edge 
e ± on which it is evaluated). 

The previous discussion may be summarized as follows. Consider dual Ising models on r,T^, which are 
finite planar or toroidal graphs; to this pair is associated a weighted bipartite graph C. For a pair of adjacent 
v e r, / e Ft, we consider X{vf)X\fv) a local observable of the pair of Ising configurations as in the left 
column of Table [l] and Y(bw) the corresponding dimer observable (in the right column). In the planar 
case, boundary conditions are as described earlier (which may affect the definition of C). By convention (for 
disorder variables, and spin flip symmetry for order variables), (Y[" = i X(vifi))r = unless the number of 
order and disorder variables is even. Let c be the local factor: 

nw e + w' e 

e£E r 



coming from normalizing the 6V weights (2.4 1. We have 



Lemma 1 (Bosonization identities). 1. In the planar case, 

I n \ I n \ In \ 

(UXivifi)) mxHfiVi)) =2c(l[Y(b i w i ) 



rt 



2. In the toroidal case, for fi, v € {0, 1}, 

V v l n \ V- v I n \ / n \ 

(lix(vifi)) (n xt (/^)) =2c(q(e A ,s B )s» A e»l[Y(b i w i )\ 

\i=l I r \*=1 / Tt \ i=l / dimer 
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Recall that in order to get consistent signs on both sides of these identities, we start from edge-disjoint 
order and disorder lines on T,T^ and deform to obtain corresponding lines on as in Figure pi In some 
cases, in order to get local dimer observables, we replace two (or an even number of) 6V disorder lines 
which almost constitute a cycle by two (or the same number of) "short" disorder lines that cross only one 
edge in M. Since the "almost cycle" encloses an even number of Ising order and disorder variables, the 
rearrangement of disorder lines leaves the sign of the correlator unchanged. 

This may be used in conjunction with Kramers- Wannier duality. We recall: 
Lemma 2 (Kramers- Wannier duality). 1. In the planar case, 

I 2n 2m \ / 1 i \ / 2n 2m \ 

(n^n^)) =2 |r| - 1 n ^ n^^n^) 

\i=l j=l I r \ee_Er / \i=l j=l / r t 

2. In the toroidal case, if [i,v € {0, 1}, 

W I 2n 2m \ / \ '*"/ 2n 2m \ 

\i=l j=l I r VeS-Er / \i=l j=l / r t 

Consequently: 

Lemma 3 (Bosonization identities: squares). 1. In the planar case, 

)2 

/ , , \ / 2ra 2m \ 

= 2 ,rt| Y[ c(ncos(^)/2)n sin (^)/2)) 

\eeEr ) \i=l 3 = 1 I dimer 

2. In the toroidal case, for fi, v £ {0, 1}, 

I 2n 2m \ \ 2 / 1 ' \ / 2 ™ 2m 

n^)IlM/i)) =-9(M^)2 |rt| II )cUs A ,e B )e^ B ]Jcos(^v l )/2)]lsinm j )/2) 

\i=l j=l / r / \ee_E r / \ i=l 3=1 

The local prefactors are not important (as long as we ensure they do not depend on the spin structure in 
the toroidal case); indeed, applying the result to empty correlators (say in the planar case) yields the (more 
appealing) : 

(2n 2m \ / 2n 2m 

Y[vM II = E dtoer JJCOS^^)^) [] Sm(0(/ i )/2) 

t=i j=i y \i=i 3=i 

Other combinations may be considered, in particular the energy e(vv') — a(v)a(v') where u and v' are 
neighboring vertices on T. Then in the simply connected case, if (//') = (vv'y, 

E r (e(vv') 2 = £^((-1)^^')-^))/^) 

The variable in the RHS is a sign depending on the parity of dimers on C crossing (vv 1 ), and may be written 
(somewhat suggestively) in terms of the height function as 1— (4>(v') — 4>(v)) 2 /2ir 2 . It is a local variable in the 
sense of dimer statistics ([30]), which relates energy correlators ( [H [5J ) to flat dimer pattern correlations 
(i). 



dimer 
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2.5 The case of the square lattice 



In the (most classical) case where T is modeled on the square lattice (ie is a portion of the square lattice, or the 
square lattice with periodic boundary conditions), the graph C is the square-octogon graph. Dimers on the 
square-octogon graph are known to map to dimers on a square lattice by the "urban renewal" transformation 
( 35J). Because of the particular importance of the square lattice and to avoid adding a mapping to the chain, 
we describe directly the mapping square lattice 6V—)- square lattice dimer (eg |21j). 

Edges of M are partitioned in, say, horizontal and vertical edges; their midpoints are the vertices of a 
square lattice C, with corresponding bipartite coloring. The position of the Ising square lattice T relatively 



to the dimer square lattice is as in Figure 10 The local rule is that if we orient a dimer from black to white 



(resp. white to black) endpoint, it makes an angle ±j (resp. ±4r) with the edge through b (resp. w), with 
6V orientation. See Figure[8] Again there is a local ambiguity around type 5 vertices, and the correspondence 



Q 



,Q 



Q 

'0' 



Figure 8: square lattice 6V configuration; associated dimer configuration 

is measure-preserving for 6V weights (wi 2 , W34, w 56 ) = (cos#,sin#, 1) and dimer weights cos(#) for NE-SW 
dimers and sin(#) for NW-SE dimers, around the generic vertex depicted in Figure [8] 

In the self-dual (critical) case (6 = 7r/4), we obtain the familiar uniform square lattice dimers. In the 
off-critical case, we have a "flipped" weighting for dimers in the sense of [9]. 

Then we have to check how electric and magnetic insertions translate in the dimer model. For electric 
insertions, the discussion is unchanged. For magnetic insertions, a small modification is needed. A sink 
defect on a horizontal 6V edge of M (with midpoint b € C) corresponds to a monomer defect at b with 
charge 1 and is still denoted by Oi(b). Similarly, 0-i(w) denotes a source defect on a vertical edge of M 
with midpoint icgCor equivalently a monomer at w. A source defect at b corresponds to a trimer defect 
on C: b is matched to two white vertices, one to its right and one to its left; this is denoted by C_i(6), as 
the total charge is —1. The weight of the trimer is the product of the weights of the two dimers it contains. 
Similarly, O\{vo) denotes a trimer at w, corresponding to a sink defect of the 6V configuration. 

With these conventions, we record the slightly modified bosonization rules in Table [2] Here v € T and 
/ G are adjacent; x £ C is the midpoint of [vf] (and may be black or white depending on the orientation 
of vf). Again, the ± signs in ±(Oi(x) — 0-i(x)) induce a minus per pair of such insertions (this matters 
only if we want to keep track of signs, which may alternatively be recovered from the asymptotic expansion). 

In the plane, monomer correlators are analyzed in [17] ; let us summarize these results in the case of 
the square lattice. If b±, . . . , b n ,wx, . . . , w n are n black and white vertices on Z 2 , the monomer correlator 
Mon(6i , . . . , b n \ Wi , . . . , w n ) is the limit of 

Z(T n \ {&!,.. ., Wn }) = (njOl(MO-lM)d"mcr 
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Table 2: Bosonization rules: square lattice 



Ising x Ising^,8V 


Dimer 


a(v)al(f) 


Oi(x) + 0_i(a;) 




±(Oi(i)-0_i(a;)) 


a(v)/J (v) 


cos(</>(«)/2) 




sin(^(/)/2) 




(O 1 (x) + O_ 1 (o:))sm(0(/)/2) 




(0 1 (x)+0_ 1 (^))cos(^)/2) 




cos(0(«)/2)sin(0(/)/2) 



as n goes to infinity, where Z(T) is the dimer partition function on the weighted graph T, and T n is a finite 
subgraph of Z 2 such that T n /* Z 2 and the uniform measure on dimers of T n converges as n — > oo to the 
maximum entropy measure on dimers on Z 2 . If b\, ... , w n have comparable pairwise distances 3> 1, we have 

Mon z2 (6 l5 . . . , «,„) ~ c m nt j|& ^ Wj|1/2 

in agreement with the heuristic interpretation of monomer defects as magnetic insertions for the free field. 
In [T7], it is also pointed out (as a by-product of the proof) that trimer defects may also be included. If we 
define 



F 

dimer 



i= l / \ /dimer 



where Si = ±1, £ i — 0> we have: 



Earner J[ (»<) ~ C m a fe J] ^ - , , 



| Ei e 3 /2 

\i=l / i<j 

where k is the number of trimers in the correlator, Vi — Vj — Q(R), R — > oo. In general, the constant a 
depends explicitly on the local geometry and the weighting conventions for dimers and trimers; with the 
the conventions in use here, a = 1. This may also be seen directly by using the 6V formulation, which is 
invariant under arrow reversal (which exchanges dimers and trimers). 



2.6 Isotropic Ashkin- Teller model 

Let us consider the case of the isotropic Ashkin- Teller model on T (see eg [5]). To each edge e of T are 
associated two coupling constants J e ,J' e - A configuration consists in a pair of spins a(v),s(v) £ {±1} for 
each v £ Vr', the weight of a configuration is 

w(a,s) = exp \p ^2 

y e=(vv')eE r J 

Notice the invariance under (s, a) O (±s, ±cr) and (s, a) ■<-» (a, s). Let us write 

where e = (vv 1 ), K e = K e (s) — J e + J' e s v s v >. Following Wegner, one may apply Kramers- Wannier duality 
to the er-field, regarding for now s as fixed. Then 

£ n^' s )= 2lrM n (e^<cosh(/^ e ) e -^) yi n emlt,tf> 

<nV r ->{±l} e e=(vv>) t:V rt ->{±1} e+ =(//') GE pt 
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where 



If we associate an 8V configuration to the spin fields (s,t) (as in Figure [2]), the 8V weights are: 



wia = e pj ' sinh(/3(J e + J' e )) w 34 = e~ pj ' cosh(/3( J e - J' e )) 

lj 56 = eP 3 ' cosh(/3( J e + J' e )) uj r8 = e~ pJ ' sinh(/3( J e - J' e )) 



Observe that: 



^56 = ^12 + W34 + W78 



Then replacing co^g, with — to^g and applying 8V duality, we obtain the 6V model with weights: 



Wl2 — W34 + W 56 + Ll>78 



sinh(2/3J e )e^ 



W12 = 



4 

—WX2 + W34 + u; 56 + w 78 



W34 = 



4 

W12 + W34 + W56 — 




W56 — 



4 



cosh(2/3J e )e /3J = 



The preceding discussion may be found in [2]. Now this argument for partition functions may be repeated 
in the presence of insertions, yielding in particular the following bosonization identity: 



for the associated 6V model with weights (cDj.2, W34, 0)56). Remark that this 6V model maps to a dimer model 
exactly when J' e = 0, ie when the Askhin- Teller model decouples into two independent Ising models. 

3 Asymptotic analysis 

3.1 Spin boundary conditions for dimers 

In order to simplify the treatment of boundary conditions, we consider the square lattice case. Let us start 
with r a simply connected portion of the square lattice, such that each face has either zero or two consecutive 
edges on the outer boundary, which is fully wired. The corresponding graph M has two types of edges: the 
regular bulk edges, and edges along the boundary. In the mapping to C, it is rather convenient to draw three 
vertices on such edges, one per quarter turn. See Figure [9j If we give weight 1 to the new edges (a pair of 
edges for each salient corner), the 6V-dimer correspondence is weight-preserving. In the isotropic case, all 
other edges have the same weight \2/2. It is then convenient to give the same weight to these corner edges; 
this has the effect of multiplying the weights of all dimer configurations by a constant (an inessential gauge 
change). 

The graph C is bipartite; we denote here Cb,Cw its black and white vertices. We are concerned with 
the dimer model on C with the special type of boundary In the analysis of the dimer model ([33 ), a crucial 
tool is the Kasteleyn operator and its inverting kernel. The Kasteleyn operator K : R Cb — > R Cw is such 
that K(w, b) = ±u>(bw) if b ~ w (uj(bw) is the edge weig ht, y/2/2 in our normalization) and K(w,b) = 
otherwise. The sign of K(u>, b) is 1 if (bw) is oriented from w to b in a fixed Kasteleyn orientation of C and 
— 1 otherwise. A Kasteleyn orientation is such that around each face, there is an odd number of clockwise 
oriented boundary edges. For definiteness, we fix the orientation of Figure [10} repeated periodically. At this 
point let us remark that in the type of boundary considered here, the (half-decimated) height is constant 
on the boundary (this is easy to see from the 6V representation); by constrast, in the somewhat easier to 
handle Temperleyan boundary conditions, the (half-decimated) height on the boundary is proportional to the 
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Figure 9: Left: graph T with a wired boundary (bold), dual graph (dashed). Right: corresponding 4-regular 
graph M. Bottom: dimer graph C. 
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Figure 10: A Kasteleyn orientation of the square lattice (dashed: a face of T) 



winding of the said boundary. However we retain a convenient feature of the Temperleyan case: a random 
walk representation of K _1 . 

Indeed, the black vertices may be partitioned into two subtypes, which are themselves the vertices of 
two dual square lattice with twice the mesh; these two types B ,Bi are represented by straight and oblique 
black squares in Figures [TUJ [TT] 

Fixing wo G Cyy, we consider h = K _1 (.,uio) as a function on Cb - (By the mapping from dual 
Ising configurations, we know that the partition function of the dimer model is positive, and thus K is 
invertible; alternatively, since the boundary height is essentially constant, Thurston's tilability criterion 
trivially applies). We will express it in terms of the Green kernel for a random walk on Bq with suitable 
boundary behavior. 

In the bulk, h is discrete harmonic except at 6o,6g, the two points of Bq adjacent to wq. Indeed, 
E b ^ w Hw,b)h(b) = S Wo by definition and Z w ~ b K(w,b)(Kf)(w) = A B J(b) in the bulk for any / G R c * 
and b £ Bo, where Ag is the positive Laplacian with edge conductances \ on the sublattice Bo- 

This may be adapted near the four types (pointing NE,NW,SW,SE when oriented counterclockwise) of 
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boundary components we are considering. In all cases, a linear combination of (K/)(u>) for w near &o (but 
not necessarily a neighbor) eliminates the dependence on B\ . The NW and SE pointing boundary arcs also 
appear in the Temperleyan case and correspond to Dirichlet (the RW is absorbed on the boundary) and 
Neumann (normal reflection). 



On a SW boundary (Figure 11 panel), by considering (Kh)(w) for w one the the circled white vertices, 
one sees that h is harmonic for the RW with transitions as in Figure [TTJ right panel. The average jump 
points in the south direction. On a NE boundary, the same argument gives that h is harmonic for a RW 
with average horizontal jump on the boundary. 




Figure 11: Left: Oblique reflection near a SW boundary component (dashed: boundary of F). Right: jumps 
of the associated RW on B . 



Lastly, we have to consider corners. Near a salient corner, h\B and h\g 1 are harmonic wrt to an 
appropriately reflected random walk on Bq, B\. However, near a reentrant corner, h\B and h\B 1 each 
fail to be harmonic at one point. Consequently, h\g Q (resp. h^) can be represented in terms of the Green 
kernel for a RW with the appropriate reflection/absorption properties on the various boundary components, 
with a pair of singularities around Wq and one additional singularity at the apex of each reentrant corner. 

In the scaling limit, one considers a fixed polygon P with boundary directions in NE,NW,SW,SE, and 
consider approximations of P by a graph T with small mesh S and the associated graph C. If we consider 
the RW on Bo with is simple in the bulk and has the transition probabilities we just described on the 
boundary, we get a sequence of processes which converges by standard arguments to a Brownian motion with 
Dirichlet /Neumann/oblique boundary condition depending on the boundary segment. Indeed, tightness is 
easily established and then one invokes well-posedness of the (sub)martingale problem for BM with oblique 
reflection in wedges ([IE]). In the cases we consider, corners are polar (as the wedge angle is at least tt/2 
and the inward oblique reflection angle is at most 7r/4, see [48]). See eg [15] for a similar reflected RW 
convergence argument. 

It is useful to estimate the probability to exit at the apex of a reentrant corner. In the continuous limit, for 
an infinite wedge, the only (up to constant) nonnegative harmonic function satisfying these oblique/Neumann 
or oblique/Dirichlet boundary conditions can be written as 5ft(z -1 / 6 ), up to rotation and centering of the 
wedge. This gives estimates on exit probabilities. For instance, in the oblique/Neumann case, this extreme 
harmonic function is 6 1 ( 1 7^ | x / 6 ) on C(0, r); it follows easily that the probability to exit the annulus A(\~ 1 r, Ar) 
on C(0, A _1 r) starting from z with \z\ = r is comparable to A -1 / 6 for A> 1. By convergence of the RW to 
the reflected BM, and reasoning on the crossings of C(0, A™ro), n E N, one obtains that the probability that 



21 



the RW started from z exits at the apex before reaching C(0,r) is > (|z|/r) _1 / 6_£ and < (|z|/r) _1 / 6+e for 
any fixed e > and \z\/r large enough. In particular, if h is discrete harmonic except at the apex and of 
order 1 on the annulus A(r, 2r), it is 0((|z|/r) -1 / 6 ~ £ ) in D(0,r). 

We want to estimate Kj 1 (.,wg), where wg — > w as the mesh S \ 0, with the direction of the edge of 
Bq containing w fixed. The standard argument is to establish existence and uniqueness of subsequential 
limits (in the proper scale). The restriction hg of Kj (.,Wg) to B is discrete harmonic except at ba,b' , the 
two vertices adjacent to wq and at one point at each reentrant corner. Fix a small e > and normalize 
hg so that \\hs\\oo,D(bi,2e)\D(bi,e) is °f order 1, where the bi's are the possible singularities. From the 
previous argument, we have hg(z) = 0(\z — c| _1 ^ 6_£ ) near a reentrant corner c, hg(z) — 0(l/\z — w\) near 
w, and hg = 0(\) elsewhere (including near salient corners), uniformly in 5. By Harnack estimates, we get 
a bound on the Lipschitz norm of hg on compact subsets of P. Along a subsequential limit, hg converges to 
a harmonic function h in P \ {w}. Moreover, h extends continuously to the boundary arcs except possibly 
at corners (this may be seen from instance by the convergence of the RW to reflected BM). 

In order to identify the limit, it is convenient (as in [39] ) to simultaneously consider hg, the restriction 
of Kj 1 (.,wg) to B>\. Up to extracting again a subsequence, we may also assume that hg converges to a 
harmonic function h on P \ {w}, which is conjugate to h. Consequently, we may write h — ip + ip, with ip 
holomorphic, and h = iip — iip. This implies that ip + ip = on a NW segment; ip — ip = on a SE segment; 
ip + ip = —vij) + it/j on a NE segment; ijj + $ = —iip + iip on a SW segment. The first two express Dirichlct 
conditions for h, h respectively; the last two are obtained by examining discrete holomorphicity at a white 
corner vertex. 

Let 4> '■ P H = {z : 5(z) > 0} be a conformal equivalence; there is a single-valued determination of 
y/4>' (as <p' is non- vanishing in P simply connected) and the RH boundary condition can be summarized as: 
ip || e~ i7TS (<p') 1/2 on the boundary. Mor eover ip(z) — 0(\z — w\ 1 ) near w, ip is of order 1 near a salient 
corner and ip(z) = OQz — c| _1 / 6_e ) near a reentrant corner c. This essentially determines ip. 

Indeed, if tp = e~ %7T ^ s y/^ip, then tpo <p~ x is meromorphic in H with at most a simple pole at (p(wo)', is real 
on K; has singularities of order 0(z~ 1+£ ) at the images of corners, which are consequently removable; and 
vanishes at infinity (since 0'(0 _1 (oo)) = 0, if we choose <p so that <p~ 1 (oo) is an interior point of a boundary 
segment). Consequently, 



Mz) = e -wy( z )i/2 ^ 



(z) - (p{w) cp(z) - (p(w) 



where ae~™/ s is the residue of ip at wq. Remark that the RHS does not depend on the choice of <f>. The 
value of the residue may be identified by substracting the inverting kernel for the full plane. Thus there is 
a unique possible subsequential limit, which gives convergence as 5 \ (for fixed w, uniformly in compact 
subsets of P \ {w}). Since we have convergence of Kj 1 (.,wg) for any sequence wg converging in P (where 
wg is always on horizontal or always on vertical edges of Bq), this also gives uniform convergence in both 
variables (w, b) when w, b are away from the boundary and from each other. 

Set e %viyV) = 1 (resp i) lib £ Bq (resp. B\) and = 1 (resp. i) if w is on a horizontal (resp. vertical) 

edge of B a . Then we have 

V \(p{z)-(p(w) <p(z)-<p(w) J J 

uniformly in (b, w) in a compact subset of P x P \ Ap. 

The reader will have noticed close similarities with expressions found in |47[ 126] ; this is expanded on 
below. 
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3.2 Monomers and Ising correlators 

We explain how to use the analysis of monomer correlators for the dimer model carried out in [17] to evaluate 
the asymptotics of Ising correlators in the planar case (for finite domains, see [51 119] ). 

Let X\ , . . . , X2m , Ui ) • ■ • ) 2/2n be fixed distinct points in C and F = <5Z 2 . We want to estimate 

(2m 2n 
i=l j=l 

as 8 \ where E is the expectation under the (unique) Gibbs measure for the Ising model on the square 
grid at the critical point, and xf — Xi — O(S) (resp. yf — Vi = 0(6)). 

For technical simplicity, let us fix a box of size n (say with sides at ±^ angles), n large enough. Let 
T" be the intersection of <5Z 2 with this box, say with wired boundary conditions. Let xf = xf + 5^^, and 
similarly yf — yf + S^ 1 ; uf (resp. vf) is the midpoint of [xfxf] (resp. [irfif]). By bosonization, we have 

(2m 2n \ 1 2m 2n \ / 2m 2n 

n n m(i^) E r t n n = (^r^w n + n - °-^)) 
i=l j=l J \i=l j=l ) \t=l j = l 

As n — > oo, the Ising measures on T (with wired boundary) and (with free boundary) converge to the 
unique Gibbs measure for Ising on Z 2 at critical temperature. Meanwhile, the dimer measure converges to 
the maximal entropy measure for dimers on Z 2 (it follows easily from the previous section). By expanding 
the RHS before sending n to infinity (and thus discarding contributions with non vanishing total charge), 
we obtain 

(2m 2n \ 2 / 2m 2n 

n^n^i)) =(-!)" e E iimcr (n^K 5 )n^(4) 

and consequently (|17j): 



j — l J &i ,fj =±i \ i— 1 J — 1 



Theorem 4. For the Ising model on Z 2 at criticality, there is c — c rn , n 7^ s.t. if the pairwise distances 
between the Xi, y } - 's are comparable and ^> 1, we have: 

(2m 2n \ 2 

Ipfa) II - E II" II to - ^r^' /2 II \y 3 - ^l 1 ^' 72 II \*i - y 3 r y] ' 2 
i=l 3 = 1 I ^,^=±1 j i^i' j^j' i,j 

Plainly, this also holds on the rectangular lattice with Z-invariant weights (see below). Let us also remark 
that, while the (pair) monomer correlation problem on Z 2 ([331 US]) has long been understood as closely 
analogous to the (pair) spin correlation problem, it is actually combinatorially equivalent to it. 

4 Relations with other approaches 
4.1 Dimer representation 

It is well-known that an Ising configuration on T can be mapped (via its high-temperature expansion) to a 
dimcr configuration on a decoration of L. It is already apparent from the discussion of spin structures that 
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there should be a direct connection between bosonization (where a pair of dual Ising configurations is in 
abelian duality with a bipartite dimer model) and this classical '"PfafHan" approach (where a single Ising 
configuration is mapped to with a non-bipartite dimer model). We now explain this connection. 

There is some flexibility in the choice of decoration: the terminal lattice introduced by Kasteleyn [25] is 
the sparsest but is non-planar; in [52], a planar decoration is introduced; a cyclically (around each vertex) 
symmetric version of the Fisher decoration is used in (3]. For our purposes, it is rather convenient to use 
(yet) another variant of the dimer mapping. 

Starting from a (planar) graph T, we consider a decorated graph Tp obtained as follows. Each vertex 
v € r of degree d is replaced with 2d vertices; one "terminal" vertex for each outgoing edge, and one vertex 
for each corner of adjacent faces. Two terminal vertices are connected (by a "road") in Tp if they correspond 
to the same edge in T; two corner vertices are connected in IV if they correspond to the same vertex v in T 
and to a pair of consecutive faces around v; and a terminal vertex and a corner vertex are connected if the 
outgoing edge corresponding to the former bounds the face corner corresponding to the latter. All in all, 
the city replacing the vertex v GT has 2d vertices and 3d edges (streets). If T is weighted, road edges of Tp 
inherit the weights and city streets have weight 1. See Figure [12} 




Figure 12: A degree 4 vertex in T; corresponding decoration in F^ 

Let us consider a polygon, viz. an even degree subgraph P of T. The corresponding dimer configurations 
m on Tp are those such that each road edge is present in m iff the corresponding edge e € T is present in 
P. It is easy to check that for any given polygon P, there are exactly 2l r l corresponding matchings of T^. 
More precisely, if m(P) is the partial matching consisting of the road edges covering an edge of P, there are 
two ways to complete it to a perfect matching within each city. 

A Kasteleyn orientation of a planar graph is an orientation such that around each face there is an odd 
number of clockwise oriented edges. The Kasteleyn operator of an oriented weighted graph M. Vr — > M. Vr is 
described by its matrix elements: K(v,v') = if v' ^ v; if v ~ v' ', K (v,v') = u(vv') if {vv') is oriented from 
v to v' and K(v,v r ) = ~uj(w') otherwise, where uj(vv') is the weight of the (unorientcd) edge (vv'). 

We want to compare the Kasteleyn operators on C and Tp, which are two weighted graphs constructed 
from the weighted graph T. For definiteness, we have to describe "compatible" Kasteleyn orientations for 
these two graphs. Let us start with a "geometric" Kasteleyn orientation for C, which is defined in terms of 
its embedding (this will also be convenient in the isoradial setting). 

For each pair v g T, / g adjacent, choose v = v(v, f) s.t. arg(/ — v) = 2v mod 2n (where we 
somewhat abusively identify the vertex v and its image in the complex plane under the graph embedding; 
equivalently, choose a square root of / — v in C) . Edges in C are either "roads" or "streets" ; let us orient 
all roads counterclockwise around vertices of T. A city vertex corresponds to a pair v £ T, f E T^; a city 
street has endpoints of type (vf), (v' f) or (vf), («/'). In the first case, the street s is parallel to (vv') and 
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is oriented in such a way that (e i M vf > +v W> , s) is direct ( viz. the oriented angle is in (0, tt)); in the second 
case, we also ask that (e % ^ v ^' +v ^ ^,s) be direct. Given the choices of i/s, this defines an orientation of 
C. Replacing viyf) with v(yf) + tt has the effect of reverting the orientation of two edges in four faces 
of C (two city faces, and the "large" faces corresponding to v and /). Thus we simply need to check the 
Kasteleyn condition for each face of C for an arbitrary local choice of v's. This is checked directly for cities. 
For a large face around v in T, let us enumerate the faces fi, • ■ ■ , fd around v in counterclockwise order and 
choose the v{vfiYs s.t. v{vfi + \) — v{vfi)) £ (0, \ ) for i = 1, . . . , d— 1 (this forces v{vfi) — v(vjd) £ (tt, )). 
Then around v we have d counterclockwise roads, d — 1 ccwise streets and one clockwise street. For a large 
face around / 6 T*, we enumerate the vertices Vx, ■ ■ ■ , Vd around / in counterclockwise order and choose the 
v{vif)'s s.t. v{v i+1 f) - v{vj)) £ (0, § ) for i = 1, . . . , d - 1 (this forces v{vf x ) - v(vf d ) £ (tt, ?f)). Then 
around / there are d clockwise roads, d — 1 clockwise streets and one counterclockwise street. We refer to 
such a choice of z/s as standard. 

We want to construct a Kasteleyn orientation for Tp derived from the orientation of C. This may be 
done as follows. Each city of Tp is centered at a vertex v ET. If v has degree d (in T), the city v around v 
has 2d outer streets and d inner streets; the face v of C around v has degree 2d. There is a correspondence 
between outer streets of v and edges of C around v; for definiteness, let us fix an edge e £ Er starting at 
v; we associate the outer street of v starting at the terminal vertex corresponding to E and going clockwise 
around v with the edge of C around v that crosses e; then the other 2d — 1 edges in C and T p are paired 
cyclically. 

In this correspondence, outer streets in Tp inherit an orientation from edges in C. Then we can orient 
each inner street in such a way that the Kasteleyn condition is satisfied for the outer triangle of the city 
containing this inner street. Since the Kasteleyn condition is satisfied for the large face of C containing v, 
it is easy to see that the inner face of the city v in Tp is clockwise odd. We still have to set the orientation 
of roads in Tp and check the Kasteleyn condition for faces of Tp which correspond to faces of T. Let us fix 
such a face / £ T^ and e <= Ep an edge on its boundary (identified with a road in IV). There is a city e 
of C associated to e. We may assume that e (oriented ccwise around /) goes from west to east; then the 
edges of e are labelled W,S,E,N. If neither or both of W and S are oriented clockwise around e, we orient 
e £ Ey f ccwise around /; conversely, if exactly one of W and S is clockwise oriented, we orient e clockwise 
around /. It is then immediate to check that if /' is the other face of T with e on its boundary, we get the 
same orientation for e reasoning around /' (since the face e is clockwise odd in C); and that the Kasteleyn 
condition is satisfied at the face of Tp corresponding to / (this uses the Kasteleyn condition for the large 
face of C corresponding to / and for the cities of C corresponding to edges on the boundary of /). The local 
situation is illustrated in Figure [13] 




Figure 13: Compatible standard orientations for C (dashed: T) and r> around an edge of T 

We want to relate the Kasteleyn operators K c : R v ° -> R Vc and K F : R Vf f -> . These are defined 
in terms of an orientation, hence in terms of a choice of v(vf)'s; changing one the ^'s results in conjugating 
Kc (resp. Kp) with a diagonal matrix with ±1 diagonal entries (in both cases, with exactly two —1 diagonal 
entries). Hence in local computations we may fix a "standard" local choice of v's as above. 

Let us start from Kp; Tp is not bipartite but its vertices are naturally partitioned evenly in two types: 
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inner and outer vertices (seen from a vertex of T). An inner vertex corresponds to a pair v £ T, f £ T^; an 
outer ("terminal") vertex corresponds to an endpoint of an edge e £ T. Let us denote by A the set of inner 



vertices and B the set of outer vertices, so that Vr F = AUB and \A\ = \B\. Naturally 



; it is 



also convenient to consider another decomposition of this space into, on the one hand, functions supported 
on A and on the other hand functions / which are "holomorphic" on B, in the sense that {Kpf)\ B — 0. Let 
us denote H B = Kp 1 (R A ) C M Vr > this subspace. 

Let us observe that Kp\^ (restricted to A and corestricted to B) has a diagonal block decomposition 
where each diagonal block (corresponding to a city) is of type 



/l -1 





Vi 









-1 

1/ 



which has determinant 2 (irrespective of the parity of its dimension d, which is the degree of the corresponding 
vertex in T). (This is for a cyclic ordering of A and B vertices around v, and a standard orientation where 
all but one of outer streets are ccwise oriented). 

Consequently, a function / £ Hb is completely determined by its restriction to B, and more precisely 
any function supported on B can be extended uniquely to a function in Hb- In particular, if b £ B, we let 
fb be the unique element in Hb s.t. (ft)\B — $b- Plainly, the /{,'s constitute a basis of Hb, and the change 
of basis from the standard basis of R Vr F to the basis {{fb)beB, {S a )aeA) (corresponding to the decomposition 
~ H B © R A ) has determinant ±1. 



Hence, if we write Kp as an operator Hb 



(viz. using the basis ((fb)beB, (S a ) aeA ) on 



the LHS and the standard basis on the RHS), we get a block decomposition 



{Kf)\1 



{Kp)\] 



Since {Kp)\ A is decomposed into (explicitly) invertible local blocks, the problem of evaluating the determi- 



nant and the inverse of Kp is essentially equivalent to the same problem for Kp d = (Kp)\^ 



Recall that 



for block triangular matrices, 



M 




N 







-M^NP- 1 
P- 1 



Let us evaluate the matrix elements of Kp. For definiteness, we work in the neighborhood of the horizontal 



edge e of F depicted in Figure 13 (with these orientations); let b, b' be the left and right endpoints of the 



corresponding road in Tp; let w be the weight of e. The values of 2w 1 fb are 

0^0 1 <- 1 



s \ 
o^o 1 -> 1 

(where we represent b,b' and the A vertices of the two cities containing b and b'; all other values are 0); 
similarly, the values of — 2w~ l fb' are: 



-2W 





t 




2G 



The extension to the general case (where the two vertices of T have general degree) is straightforward. We 
see that Kpfb and Kpfv are supported on the four A-vertices adjacent to b or b'; let us label this vertices 

\ s 

I b^b' t 

S \ 
0,2 a 3 

Then the restriction of Kp to columns corresponding to b, b' and rows corresponding to a\, . . . , 0,4 reads: 

1 1 —w w x 
w —w 1 1 

Using the trigonometric parametrization w = tan 0/2, we notice the identity: 



/cos(0/2) -sin(0/2)\ / 1 1 -tan(0/2) tan(0/2)\ _ , , . /cos(0) 1 

\sm(e/2) cos(0/2) ) ^tan(0/2) -tan(0/2) 1 1 J ~ C ° S( 1 ' \sm{9) 



sin(6») 
cos(0) 1 



If we consider T F a Fisher decoration of 1^ (with dual weights (2.1 1), there is natural 1-1 correspondence 



between type A vertices of Tp, T F and a natural 2-2 correspondence between type B vertices. At this stage 
it is easy to see that K F may be written in terms of Kp, yielding a version of Kramers- Wannier duality at 
the Fisher representation level. 

Let us turn to K C - As C is bipartite, R Vc ~ R v o © R v c where V§ (resp. V^) designates the black 

/ K BW \ 

(resp. white) vertices of C . The corresponding block decomposition of Kq reads [ ^wb q j with 

Kg w : R v c _> M,Vc r and K™ B = -{K^ B Y, Pf(K c ) = det(K BW ) (where Pf designates the Pfaffian of the 
antisymmetric matrix Kc in a standard basis); this is the standard for planar bipartite graphs. 

Here, we can establish a natural bijection between black vertices of C and type A vertices of Tp; and 
between white vertices of C and type B vertices of Tp-. This gives the following labels for vertices of C near 
the edge e: 

b" a 4 

\ / 
ai — ► b' 

I t 
b -> a 3 

a 2 b" 1 
Then the submatrix of K BW with columns indexed by a\, . . . , and rows indexed by b, b' is 

cos(0) 1 -sin(0) 0\ 

sm(6) cos(6») \) [ W) 

We conclude that there is a rotation matrix O (block diagonal with 2x2 blocks) and D diagonal such that: 

OKp = D ■ Kg w 



where the entries of D are of type cos(0 e /2), with 8 e such that w e — tan(# e /2). Remark that the Pfaffian of 
Kp enumerates (weighted) dimers on Tp (corresponding to polygons on V), and the determinant of K BW 
enumerates (weighted) dimers on C. 

Let us point out some consequences. If the graph T is periodic (viz. there is a group of isomorphisms 
~ Z 2 operating on T with a finite quotient, say generated by a right shift R and an up shift U), Tp and C are 



27 



also periodic and we may choose orientations periodically. Plainly, Kp : Hb — > R A and Kq W : M. v c — > M 17 *?' 
and the connecting matrices O, D commutes with translations. If z, w € C*, and V is a (complexified) vector 
space on which R, U operate, let V Z:W = {v € V : Rv = zv, Uv = tow}. Then (Kp) z _ w ■ (Hb)z,w — > (C A ) ZjU) 
and {K% w ) z , w : {C v c) z , w -> (C^)^ are related by O p (£ F )* )U; = # P • {K% w ) z ^, where O p ,£ p do not 
depend on z,w. Consequently, if we consider the characteristic polynomials Pp(z,w) = det((Kp) ZtW ) and 
Pc(z,w) — det((K c ) z ,w) (this is well-defined modulo sign for a choice of standard bases for the various 
vector spaces), we have 

P F {z,w) = cP c (z,w) 

where c ^ depends only on the degree distribution and the weights (the block (Kp)\^, operates only within 
cities and contributes ±2l rp l to the constant c, where T p is a fundamental domain of T). See [3^1 [HI 130] for 
related questions. 

A setting where T is not necessarily periodic but is still tractable is the isoradial case. Start from a tiling 
of the plane by rhombi (the diamond graph <^>), say with edge length 1 (for normalization). This is a bipartite 
graph. By retaining every other vertex and connecting them by a diagonal in each rhombus, one obtains an 
isoradial graph V (its dual is also isoradial and corresponds to the other half of the rhombus tiling vertices) . 
See Figure [14] There is a special choice of Ising edge weights on T associated to its isoradial embedding: the 




Figure 14: a rhombus tiling <) (thin); associated pair of dual isoradial graphs T, (bold, dashed) 

critical Z -invariant weights (^HlEl ?]), which are in particular significant from the Yang-Baxter point of 
view. If e € E-p, its weight is 

w e = e~ 2J ' = tan(0 e /2) 

where 9 e is half of the angle under which e is seen from the center of a face of T with e on its boundary. 
Remark that on the dual graph, u> e t = tan(? — ^), a parameterization of the Kramers- Wannier duality 



relation (2.1 1 



Let us consider C, which is up to now embedded rather arbitrarily. Each pair of a black 6 and white 
w vertex connected by a road in C corresponds to a pair v € T, f E T', ic to an edge in <); we choose the 



coloring of C s.t. (vf,bw) is direct (as in Figure 13). Let us embed both b and w at the midpoint of this 



edge of <) (this is not a proper graph, but is nonetheless convenient here). Then C is itself isoradial, and 
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corresponds to the diamond graph <)' obtained from <£> by dividing each face of (} in four isometric rhombi, 
and inserting a fiat rhombus on each edge of <). This is somewhat degenerate but preserves most of the 
isoradial machinery (in particular the all- important discrete exponential functions, |42l [32j [7] ) . See Figure 
US] 




Figure 15: a rhombus tiling ^>; associated rhombus tiling <0>' and bipartite isoradial C (with flat rhombi 
blown up for legibility; to be glued back along dashed lines) 

Associated to any bipartite isoradial graph Co (with diamond graph <%) is a Dirac operator K introduced 
in |32j . It is a complex weighted adjacency matrix with matrix elements: K(w, b) — i(u — v) if w, v, b, f are 
vertices on a face of <Oo listed in counterclockwise order, with b £ Vq , w € Vq (adjacent black and white 

vertices). This defines K : C c « ->• C c ° . 

In the special case of <0>' and C, for each pair v G T.f e T*, we have chosen an angle v{vf) s.t. 
/ — v = |/ — v\e 2w ^°^' . The pair (v, f) corresponds to a pair b, w of vertices in C; set v{b) — v(w) = v{vf). 
Then: 

K(w,b) = ie l (^ +u(b »Kc(w,b) 

where Kc is the (real) Kasteleyn operator we have been using and K the critical Dirac operator. This is a 
simple direct check (distinguishing the cases of roads and streets) , which informed the choice of Kasteleyn 
orientation on C . 

4.2 Fermionic variables 

For the reader's convenience, we include a discussion of fermionic variables and associated correlators in 
some of their various incarnations. 

Following Kadanoff and Ceva j^Tji we have considered a fermionic variable as the product of adjacent 
order and disorder variables: ip(vf) — a(v)/j,(f) (which is invariant by Kramers- Wannier duality). This leaves 
an ambiguity on the sign of a correlator containing tp, which is insufficient for some local computations. 
This ambiguity may be resolved in a couple of different ways; let us indicate one which is immediately 
compatible with the previous discussion, restricted for simplicity to the case of a pair fermionic correlator 
(ifi(vifi)'ip(v2f2)} in a simply connected domain. 

In order to assign a sign to this correlator, it is enough to specify a disorder line from fi to /2, or rather 
an equivalence class of disorder lines. This may be done in terms of the choice of angles v(v, /). In a slightly 
more general setting, one can consider an assignment (v, f) H> 2v(v 1 f) 6 K s.t. arg(/ — v) = 2v{v, /) 
mod 2-7T. Given v(v\, f\) and a simple path 7 : v\ — > V2, one can define a natural parallel transport 
along 7 by the following local rules: there is a sequence ((v k , f k ))k=o,...n such that (v°,f°) = 
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(resp. (v n , f n ) = (v 2 ,f2)) and for all k, ((v k ,f k ), (v k+1 , f k+1 )) is of one of two types: ((v, /), (V, /)) with 
v,v' consecutive vertices around / in clockwise order; or ((«,/), (v, /')) with /, /' consecutive points on 7 
and v the vertex on the lefthand side of //' £ E T i ■ We may define sequentially angles v'(v k f k ) so that 
v'iv / ) — v(vqJq) and the increment v' '(v k+1 f k+1 ) — v'(v k f k ) is in (— 7r,0) in the first case and in (0,7r) in 
the second case (with arg(/ fc — v 1 ) = 2v'(v k f k ) mod 2-7T for all k). Let us denote T^v^ifi) — v'{v n f n ). It 
is easy to check that T^rT^v = v — ir if 7 r : / 2 — > fx is 7 taken with reverse orientation. 

Given a choice of z/s, we may specify the sign of (ip(vifi)'ip(v2f2)) by evaluating this correlator with 
respect to a disorder line 7 : /1 — > /a satisfying: 1^^2/2) = T 7 v(vxfx) mod 27r; in other words, if 
(il>(vifi)ip(y2f2))'y is the evaluation of the correlator wrt 7, we may set 

W«l/lM«2/ 2 )) = e i(t/ ' (, ' 2/2) " T ^ (, ' l/l)) (^K/i)V'(«2/2)) 7 

Remark that with these conventions, the ■(/''s are (as they should) anticommuting variables, viz. (ip(vxfx)4 ! ( v 2f2 
-M«a/2M«i/i)>. 

Let us consider an Ising model on a simply-connected 1^ (spin variables are defined on faces) in its low- 
temperature expansion (a polygon on T). If (vifi), i = 1,2 are two vertex-face pairs, (f 2 )<r^ (fx) = (—1)™ 
where n is the number of polygon edges crossed by any path 7 : fx — > /a (n mod 2 does not depend on the 
choice of 7, since the polygon has even degree at each iieT). Disorder variables (which are here located on 
Vr) have the following low-temperature representation: a polygon with disorders at V\, V2 is a subgraph of T 
with even degree at every v € Vr \ {vi, V2} and odd degree at v\, «2- Then the correlator (V 1 (/ivi)i/' 
may be represented in terms of polygons with disorders at v\ , V2 and counted with a sign depending on the 
parity of the number of polygon edges crossed by 7 : fx — > fa (which does not cross a specified disorder line 
V\ -> v 2 ). 

From the low-temperature expansion, one may move to the Fisher representation. From a polygon on T 
with disorders at v\, . . . , v 2n , one can define a matching on the Fisher graph . Each polygon edge descends 
to a matched road on ; within each regular (^ vx, ... , «2n) city there are two choices to complete the 
matching. In each disordered city (corresponding to one of v\, . . . ,v 2n )), remove an inner (A-type) vertex 
corresponding to the pair (vifi). Then it is easy to check that the road matching may be extended in a 
unique way in the city with a monomer at ai ~ (vifi). Consequently, in terms of the dimer configuration on 
Tp, the correlator ( - i/ , (^i/i)V'( t; 2/2)) is expressed (up to a multiplicative factor 2 2 coming from the disordered 
decorations) as the partition function of matchings on Tp \ {vi, . . . , V2 n } counted with a sign depending on 
the parity of matched edges crossing a path /1 — > /2. As is well-known, this is (up to sign) the inverting 
kernel of Kp evaluated at the pair ai, a 2 : Kp 1 (ai, 02). 

As discussed earlier, Kp 1 may be expressed in terms of K^, 1 by linear algebra manipulations. Alterna- 
tively, by bosonization the correlator (?/'( u i/i)V'(' l '2/2)) may be expressed in terms of magnetic/electric dimer 
correlators on Kc, which in turn are readily identified in terms of the inverting kernel K q . 

A variant of the (low temperature) polygon representation uses edge defects: every vertex has even 
degree, but in some prescribed edges one half-edge is present in the polygon and one half-edge is absent. 
Consider the partition function of these polygons with defects at e\ = (viv[),e2 = (v2v' 2 ) € £r, where as 
before configuration are counted with a sign depending on parity; this is essentially the spin Ising observable 
used in [5]. By erasing the defective half-edges, one obtains a polygon with one vertex defect at V\ or v[ and 
one vertex defect at V2 or v' 2 , where the edges ei,e2 are vacant. It is then easy to see that this observable 
may be expressed as a linear combination of pair fermionic correlators. 

As explained in [43], the correlator (V , (^i/i)V'( u 2/2)) has a natural representation in terms of the random- 
cluster representation of the model (or rather the associated fully-packed loop representation), in the simply 
connected, wired boundary case. In the case of a simply connected domain with a wired and a free boundary 
arc, if (wi/i) is taken "in the bulk" and (^2/2) is at one of the boundary condition change point on the 
boundary, one obtains the FK Ising observable of [H>] (up to normalization). 

An important property of fermionic correlators {4>(vf)X) is that they satisfy local linear relations. In |14) . 
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it is shown (on the square lattice, see [32] on isoradial lattices) how to obtain a relation on the (-0(v/i)X)'s, 
where the fi's are the faces adjacent to v (and by duality on the (4>(vjf) ■ ■ •), where the v/s are the corner 
of /). Here X represents arbitrary (fixed) order and disorder insertions away from vf. Linear relations of 
this type are instrumental in determining the scaling limit of FK and spin interfaces |46j . 

These relations may be derived in any of the various representations of fermionic correlators; let us sketch 
one possible way for the Kadanoff-Ceva representation. Consider: 

F(vf) = (^(vf)X) 

where X represents order and disorder insertions away from vf (say more than one lattice spacing away; 
order variables are on faces in order to simplify the comparison with the Fisher representation). Let 
(vv') be an edge of T and (//') be the dual (oriented) edge in T^; we want to relate the four values 
F(vf), F(v' /), F(vf), F(v' /'). Let us consider partial partition functions Z± t ± where 



Z e , e , = 



£ 

r(f)=eMf> 



X(a) 



n 



e=(^')Ei5 rt 



Remark that if we choose disorder lines, a pair of disorder variables can be represented as a (nonlocal) 
random variable. For dcfmitcncss, assume that (vv 1 ) goes from left to right, (//') from bottom to top; we 
choose a disorder line 7 from v to v" (another disorder variable included in X) which does not contain (vv'); 
and a disorder line 7' starting from v' obtained by concatenating (v'v) and 7. We use a local choice of z/s 
as in Figure 13 (F(vf) = (vf)X) 7 , where ± depends on the choice of f's and 7). Then 



where 
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Each pair (vf) corresponds to a black vertex of C; hence we may see F as an element of M. v c , and we have 



just checked (see (4.10)) that K^ w F = (away from other insertions, which create monodromies and/or 



poles) . 

More generally, if / € 
have: 



v c and KS W f = in some region, then around an edge e G Er as in (4.101, we 



ie (f(ai 



ff(ai)\ 

cos(0) 1 -sin(0) 0\ f(a 2 ) _ [0 
sm(0) cos(9) 1 ) f(a 3 ) ~ [o 

, f((Xi)) lies in a two-dimensional real vector space, which may be identified with 



via: 



(^(z e e- lu{ai) ), . . . ,3?(z e e- w ^)) 



Here 6 € (0,tt/2) 
we have g : E-p — > 
M-linear relation: 



5 arg and (v(ai), . 



■ ,v(a 4 )) = (|,7T 



2' 2 



e 37r 

2' 2 



j). Hence associated to / 



such that if e, e' £ Ep are the two edges of T corresponding to a € Vq, we have the 



R(ze/e-*<°>) - f(a) = ^(z e e- 1 ^) 



which is (up to a rotation of the z's) the notion of S -holomorphicity ([5]). 
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4.3 Quadratic relations 



In |45j . Perk identifies quadratic identities between spin correlators, which may be interpreted as discrete time 
Toda equations. Let us briefly indicate how to recover such identities (on a general graph) from bosonization. 

We consider dual Ising models on (simply-connected) T, F^; for i = 1, 2, e,* = (viv'j) € Er and e\ = (/«/,')■ 
Let us start from 

which by Kramers- Wannier duality is 

c « < 7(« 1 )^(« a )) r <M t (i/i)/* t («2)>rt - M« l )«T(^)>r<M t K)M t («2))rt) 

where throughout c is an (explicit) product of local factors (independent of insertions). By bosonization, 
this is rewritten as: 

c ((^iM^M^iM^VKM^W - (M w i)M4)mK)m(«2)ct(i4)o-M)6v) 

The edges ei,e2 € -Er are identified with a pair of vertices on the 6V graph M. In both cases, we have a 
6V configuration with disorder lines from v\ to v[ and v-i to v 2 , which may be drawn so as to intersect two 
consecutive edges of M around ei, e 2 . A direct examination shows that terms corresponding to a sink and a 
source around e\ and a sink and a source around ei are counted with the same sign in both correlators; this 
leaves terms corresponding to two sinks around e\ and two sources around e 2 or vice-versa. Remark that 
having two sinks on edges around e\ forces the 6V type of ei; and that such a pair of edge defects may be 
alternatively represented by a vertex defect, where the arrow configuration around v± is an 8V sink. 

Hence the above combination of spin correlations may be expressed as a 6V correlator with a pair of 
(magnetic) 8V defects at ei,e 2 counted with a sign (electric) depending on the height variation from, say, 
v[ to v\. This is plainly symmetric in T H and consequently: 

(CT(ui)o-(v2))r<o-K)< 7 K))r - ( cr ( u i)°"( w 2 )>r(o'(wi)cr(w2))r = 
c(^ t (Ay(/ 2 )) r t(at(/0at(^)) rt -^t(/0 f Tt(/^)) rt (at(/{)at(/ 2 )) rt ) 

where c is an explicit local factor depending on the choice of normalization. 



4.4 Consequences for dimers 

Bosonization allows to transfer questions about Ising correlations into problems on (bipartite) dimer corre- 
lator. The latter are typically analytically easier to handle and come with a natural free field interpretation. 
Let us point out however that some Ising results may be transferred back to dimers; we shall mention a few 
of these. 

For the dimer height function on the square lattice, we know |31j that the scaling limit is the free field, 
which in the plane has the distributional invariance <p ■<-> —<j>. However at the discrete level this is not 
immediately apparent. An easy way to see it is to map the dimer configuration to a 6V configuration, 
use the arrow reversal of the latter, and map back to a dimer configuration. The mapping is (almost) 
deterministic and (almost) involutive, with ambiguities coming only from type 5-6 vertices (Figure [8]). 

In [T7j, the asymptotics of electric correlators of type (exp(ia(4>(y) — <p(x)))) are evaluated for a charge 
a G (— J, 4), the limiting case a = \ appearing as rather delicate. On the planar square lattice Z 2 (with 
critical Ising weights), we may write 

(2n \ / 2n \ / 2n 
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where u = ^p- and (Z 2 )^ ~ u + Z 2 . Using bosonization on both sides yields: 

/ 2n \ / 2n 

Edimor Y[(0! + 0-l)(Vi + |) J = Edimcr [YlcOS^V,) /2) 



\i=l / \i=l / 

where the dimer graph is |(u+Z 2 ). This is an exact identity between electric and magnetic dimer correlators, 
and is a discrete version of (a by-product of) T-duality for the limiting free field. Since the monomer 
correlators are worked out [17], this also allows (for the square lattice) to evaluate electric correlators with 
half-integer electric charge. 

Lastly, let us consider the case of a subgraph C of (5Z 2 (carrying dimers) corresponding to an Ising 
model on <5(^ + 21?) with wired (or dually free) boundary conditions along (a (^-approximation of) a fixed 
smooth closed loop 7. On the one hand, the inverse Kasteleyn matrix K^ 1 (.,.) may be identified as a 
(properly normalized) pair fermion correlator. On the other hand, such pair correlations (under a different 
representation, see above) are analyzed in [26J , building in particular on the discrete integration of squares of 
S'-holomorphic functions. This shows that the asymptotics of K^, 1 derived above (in an a rather pedestrian 
fashion) for polygonal domains actually holds for a much larger class of (simply-connected) domains. Let us 
repeat that this choice of boundary corresponds to a flat boundary height for dimers. Then convergence of 
the height field to a free field with constant (Dirichlet) boundary conditions follows easily from Sects 3-5 in 

[UJ. 
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